function data = gcpick_readdatafile(filename)
tem = readsac(filename);
data.station = [tem.KNETWK, '.', tem.KSTNM];
data.wave = tem.DATA1;
data.lat = tem.STLA;
data.lon = tem.STLO;
data.stel = tem.STEL;
data.reftime = [tem.NZYEAR, tem.NZJDAY, ...
    tem.NZHOUR, tem.NZMIN, tem.NZSEC, tem.NZMSEC];
data.axis = tem.KCMPNM(3);
data.delta = tem.DELTA;
data.npts = tem.NPTS;
data.flag = 0;
end